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Heavy Flavour Lifetimes and Lifetime Differences 

Jonas Rademacker a 

a CDF Experiment, Fermilab, and University of Oxford 

We give an overview of heavy flavour lifetime measurements, focusing on recent results from the Tevatron and the B factories. 



1 Introduction 

In the first part of this article we summarise the status and 
latest measurements of B-hadron lifetimes and lifetime ra- 
tios, including some recent result from the Tevatron and the 
B factories, and compare those results with the predictions 
from Heavy Quark Expansion (HQE). Future prospects for 
lifetime measurements at the B factories and the Tevatron 
are discussed. 

In the second part, we review the status and prospects of 
measuring the difference between the lifetimes of the two 
CP eigenstates in the B^ -B^ system. 

2 Lifetimes and Lifetime Ratios 

2.1 Theoretical Predictions on B Hadron Lifetimes 

Life time measurements in the heavy quark sector gain spe- 
cific significance due to the precise predictions of Heavy 
Quark expansion (see e.g. [1], [2]) thus providing a test- 
ing ground for this theoretical tool that is frequently used, 
for example to relate experimental measurements to CKM 
parameters like T d to \V c b\ or Am^/Am^ to \V ts /V t d\. 

The hierachy expected for b hadron lifetimes is [3]: 
r(B c ) « r(E°) 

~ r(A b ) < t(B°) ~ r(B°) < r(B") 
< r(S") < T (Ob). 



Recent HQE predictions for the lifetime ratios are [4] : 

• t(B")/t(B°) = 1.06 + 0.02 

• t(B b )/t(B°) = 1.00 + 0.01 

• T(A b )/r(B°) = 0.90 + 0.05 



2.2 The B Factories 

The B factories BaBar and Belle at the asymmetric e + e" 
colliders PEP-II and KEK have collected 123fb _1 and 
145 fb" 1 worth of data respectively up to May 2003, run- 
ning at the T(4S) resonance. 



opposite B 
a z(rms)= 1 ' 



Y(4S) 





reconstructed &Y 
az(rms)=70|H 



Figure 1. Schematic of lifetime measurements at BaBar 
(from [12] with added comments). 



2.2.1 Lifetimes at the B Factories: Method 



The T(4S) decays to B°, B° or B + , B", nearly at rest in 
the T(4S) rest frame. By colliding e + and e" of different 
energies, the CM frame is boosted (fiy ~ 0.54 for BaBar 
and ~ 0.43 for Belle) such that the B^ and B^ travel a 
measurable distance in the detector before decaying. 



Because the B mesons decay virtually at rest in the T(4S) 
frame, their momenta in the lab frame are known from 
the beam momentum. This constrains the decay dynam- 
ics considerably with the important consequence that the 
decay vertex of a B meson can be obtained from a single 
decay product, by intersecting its track with the beam axis. 
The decay distance along the beamline (z) is directly pro- 
portional to the proper decay time for a given beam mo- 
mentum (small corrections apply [12]). Lacking primary 
vertex information, it is the distance between the decay ver- 
tices of the two B mesons that is used for measuring the life 
time. This distance is typically - 250 jj,m and ~ 200 /im at 
BaBar and Belle respectively. 

In the standard method, one B meson is fully reconstructed, 
(B rec ), and another one partially, from as little as one or two 
tracks (B opp ), with a correspondingly degraded vertex reso- 
lution. The lifetime difference is calculated from the differ- 
ence in the position along the beam line (Az = ZB reco _ z pp) 
of the two B vertices. This is illustrated in Fig. 1. The res- 
olution function is modeled using Monte Carlo simulation. 
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Figure 2. Az Resolution at BaBar for B + -> J/^K + [12], with 
exp(Az/250yL/m) superimposed for illustration. 

In Fig. 2, the Monte Carlo generated resolution function for 
Az at BaBar for the decay B + — > J/^K + is shown [12]. An 
exponential with a mean decay distance of 250 mum, rep- 
resenting approximately the signal distribution before de- 
tector effects, is superimposed, illustrating how the signal 
is of a similar width as the resolution function, which must 
therefore be modeled carefully. This modeling of the res- 
olution function, together with the modeling of the back- 
ground distribution, is the biggest systematic uncertainty 
in both experiments. The so-called "outliers", a relatively 
small number of events with very large reconstructed Az, 
represent a particular problem. Both experiments are able 
to control it well enough to keep the systematic error below 
the statistical uncertainty. 

Both experiments describe the At = Az/ (c (J3y) r ) distri- 
bution in terms of three components: signal, background 
and outliers. The beam constrained mass (shown in Fig. 3 
for Belle) is used for an event-by-event signal probabil- 
ity. The fraction of outliers is a free fit parameter. The 
B^ and B + distributions are fit simultaneously in an un- 
binned likelihood fit. Besides these commonalities, there 
are some differences in the event selection and modeling 
of the resolution function which are described in detail in 
the publications by the respective experiments [11] [13]. 

2.2.2 Results 
Belle 



Using the following fully reconstructed hadronic decays: 

B o -> D ( * } -(7r + ,p + ), J/<AK°, J/^K*°, B + -> dV, J/<AK + , 
Belle find the following B^ and B + lifetimes [13]: 

r B o = 1.554 ±0.030 + 0.019 ps 




5,2 5,22 5,24 5,26 5,23 5.3 



Figure 3. N- of events vs Beam Constrained Mass at Belle [13] 

r B + = 1.695 ±0.026 + 0.015 ps 
Tb+ /t b o = 1.091 ±0.023 ±0.014 



The fit result to the data, showing separately the back- 
ground and the outlier contribution, is shown in Fig. 4. 

BaBar Fully Hadronic 

Using the following fully reconstructed hadronic 
decays: B° D ( * } -(7r + ,p + , a+), J/<AK°, J/^K* and 

B + dV,J/^K + ,i/<2S)K + BaBar find the following 
B° and B + lifetimes [11]: 

r B o = 1 .546 ± 0.032 + 0.022 ps 
t b + = 1.673 ± 0.032 ± 0.023 ps 
t b+ /t b o = 1.082 ±0.026 ±0.012 



The fit result to the data is shown in Fig. 5. 
More results from BaBar 

As mentioned earlier, the decay kinematics at BaBar and 
Belle allow to find the z position of a decay vertex from as 
little as a single track. While in the previously mentioned 
measurements, one of the pair of B's is fully reconstructed, 
BaBar also published a set of measurements where also 
the B rec is reconstructed partially. These are summarised 
in Table 1 . 
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Figure 4. Belle's life time fit [13]. 
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Figure 5. BaBar's life time fit [11]. Solid line: total fit. Single 
hatched: total background. Cross hatched: outliers. 





B u -> D*" (partial^ v [14] 




1.529 ± 0.012 ± 0.029 ps 




B u -> D*-(partial)(7r + ,p + ) [15] 


T B o 


1.533 ±0.034 + 0.038 ps 




B u -> D*~€ + v [16] 


T B o 


1.523^;^ + 0.022 ps 




Di-lepton (prelim) [17] 


T B o 


1.557 ±0.028 + 0.027 ps 


T B + 


1.655 ± 0.026 ± 0.027 ps 


T B +/T B o 


1.064 ± 0.031 ± 0.026 ps 



Table 1. BaBar's results from partially reconstructed decays. 
Here, a "partial D*" is a D* decaying to to D°7r, reconstructed 
using kinematic constraints and the pion momentum, only, with- 
out reconstructing the D° [15]. 

^t(B + )/t(B°): 
4- EXP: 1.073 ± 0.014 

■-h HQE: 1.06 ± 0.02 

T<B s )/t(B°): 

^ EXP: 0.949 ± 0.038 

i-h HQE: 1.00 ± 0.01 

T(A° b )/T(B°): 
• — i — i EXP: 0.797 ± 0.053 

~ . HQE: 0.90 ± 0.05 

0.7 0.8 0.9 1 1.1 1.2 

Figure 6. Status of lifetime ratio measurements, incl. results 
from BaBar, Belle, CDF, LEP. r(B s )/r(B d ) and r(B + )/r(B d ) 
from [18], r(A b )/r(B d ) from [19]. 

2.2.3 Status of Lifetime Measurements Including Re- 
sults from BaBar and Belle 

Since the B factories have started taking data, they have re- 
duced the error on r B + /r B o by half. Fig. 6 shows the current 
status of the life time measurements and compares them 
with HQE predictions. The r B + /r B o measurement, domi- 
nated by the precise results from the B factories, is already 
more precise than that of the HQE prediction, and we can 
expect further improvements in the near future. 

The situation is different for the B s and the A^, which are 
not accessible at the B factories. The experimental pre- 
cision of the r B o /r B o measurement lags behind that of the 
HQE calculations. For the At, experiment and theory don't 
appear to be in very good agreement, but the experimental 
and theoretical uncertainties are still rather large. Improved 
measurements and calculations are needed for clarification. 

Both, B s and A b particles are produced abundantly at 
hadron colliders, from where we expect dramatically im- 
proved lifetime measurements in the near future. The 
hadron collider currently producing large numbers of B s 
and Ab is the Tevatron at Fermilab. 
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Figure 7. Luminosity at the Tevatron 



LI: 2 XFT tracks, p t 
Pa + Pa > 5.5 GeV. 

L2: 
2-body: 
e.g. B° 7in 


> 2 GeV, A0 < 135°, 

Multi-body: 
e.g. B° ^ D s tt 


100yum< IP < 1mm 
20° < A0 < 135° 
L xy > 200 yum 
IP of B < 140 yum 

L3: Same with refined track 


120yum< IP < 1mm 

2° < A0 < 90° 
L xy > 200 yum 

s & mass cuts. 



Table 3. The CDF hadron trigger. A0 is the angle between the 
tracks in the transverse plane. IP is the impact parameter in that 
plane. L xy is the decay length in the transverse plane, which can 
be calculated from the impact parameters and A0. 



Projected / £dt / (ftr 1 ) 



Year 


Baseline 


Stretch 


2002 


0.08 


0.08 


2003 


0.2 


0.32 


2004 


0.4 


0.6 


2005 


1.0 


1.5 


2006 


1.5 


2.5 


2007 


1.5 


3.0 


2008 


1.8 


3.0 


Total 


6.5 


11. 



Silicon vertex trackers providing excellent proper time res- 
olution, sufficient to resolve the expected fast oscillations 
in the B^ system. The excellent impact parameter resolu- 
tion is used for triggering on B -events. Both experiments 
have increased their muon coverage since Run I, and have 
an efficient di-muon trigger for finding B° — > J/ifsX decays. 
D0's //-trigger covers a particularly large pseudo rapidity 
range up to \rj\ = 2. 

IP Trigger 



Table 2. Projected integrated luminosity at the Tevatron for base- 
line and best-case ("stretch") scenario. The total integrated lumi- 
nosity by 2008 is expected to be between 6 fb -1 and 1 1 fb -1 . 

2.3 Lifetimes at the Tevatron 

2.3.1 Run II 

CDF and D0 have been taking data at Tevatron Run Ha for 
about two years. For pp collisions at 1.96 TeV, the bb pro- 
duction cross section is cr bb -0.1 mb. The integrated lu- 
minosity delivered until June 2003 is shown in Fig. 7. The 
integrated luminosity at Run Ha is expected to be 2 fb -1 . 

The projected luminosity for each year until 2008 is listed 
in Table 2, for two scenarios: The base-line scenario, and 
a best-case scenario ("stretch"). The total integrated lu- 
minosity at the end of Run II in 2008 is expected to lie 
between 6 fb -1 and 1 1 fb" 1 . 

2.3.2 D0 and CDF 

Both experiments at the Tevatron have undergone major 
upgrades for Run II, optimising their B physics potential. 
The most significant upgrade at D0 is the introduction of 
a magnetic field and a new tracking system providing pre- 
cise momentum information. This significantly improves 
the mass resolution. CDF also improved its tracking with 
a new, faster drift chamber. Both experiments have new 



One of the most innovative improvements for B physics 
at the Tevatron is the large-bandwidth hadron trigger at 
CDF, which triggers on the impact parameters of tracks 
at Level 2. The eXtremely Fast Tracker (XFT) uses pat- 
tern matching to find tracks in the COT (drift chamber) 
within 5.5 yus, with about 96% efficiency for momenta 
above 1.5 GeV. These XFT tracks are combined with 
tracks in the Silicon Vertex Detector by the Silicon Vertex 
Tracker (SVT), which makes impact parameter informa- 
tion available at Level 2 to a precision of ~ 50 jim. The 
2-Track hadron trigger combines the information on the 
direction (XFT), momentum (XFT) and impact parameter 
(SVT) to trigger on hadronic B decays. The trigger re- 
quirements for the two scenarios, 2-body and multi-body B 
decays, are given in Table 3. The SVT+lepton trigger for 
semileptonic B decays has impact parameter requirements 
on one track only and requires additionally an electron or 
muon with p t > 4 GeV. 

D0 also has impact parameter information available at 
Level 2, and will have a lepton+displaced track trigger, 
which was however not yet available for the data presented 
here. 

For lifetime measurements it is essential that the bias due 
to the impact parameter cuts in the trigger is corrected for. 
We will first consider measurements that do not suffer from 
such a trigger bias, and then those that do. 
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Figure 8. Projection of Fit to cr B + from — > J / i//(pi + \i )K + at 
CDF, using 70 pb" 1 of data. 

2.3.3 Measurements Using Fully Reconstructed De- 
cays, Without IP Trigger 

Both experiments have published results from fully recon- 
structed hadronic B —> J/if/X decays from the dimuon trig- 
ger, which are not biased by any impact parameter cut. An 
example fit (CDF, 70pb _1 ,B+ -> J/if/Qi + fi-)K + ) is shown 
in Fig. 8. The signal is modeled with an exponential, the 
background by a prompt component and two positive and 
one negative exponential tails (only one positive tail for B^ 
because of lower statistics). Signal and background func- 
tion are convolved with a single Gaussian to take into ac- 
count detector effects. The B mass is fit simultaneously 
and provides an event-by-event signal probability. D0 
use a somewhat different approach, as illustrated in Fig. 9, 
modeling the background using a separate fit to the right 
sideband. The left sideband has a long-lifetime component 
from incompletely reconstructed other B decays. This B 
contamination in the signal region is modeled from Monte 
Carlo and found to be 12%. 

The results are given in Table 2.3.3. The table shows that 
the error on the life time ratios obtained from B — > J/if/X 
decays is about twice that achieved in Run I, all channels 
combined. By the end of this year, CDF is expected to have 
collected ~ 300 pb -1 , four times as much as used for the 
analyses presented here, so we can expect CDF to achieve 
the combined Run I precision using the exclusive channels 
alone by the end of this year. 

2.3.4 Measurements Using Partially Reconstructed 
Decays, Without IP Trigger 

Inclusive r B 

Since all B hadrons can decay to J /if/, reconstructing J /if/ 
vertices allows to find an average B lifetime, where the 
composition of the sample depends on the detector and se- 
lection criteria. Since the decay is not fully reconstructed, 



Absolute Lifetimes (D0 [20] and CDF [21], Run II pre- 
lim.) 



B u + 


-> J/<A0uV)K + 


D0 


1.76±0.24(stat)ps 


B u + 


-> WOu V)K + 


CDF 


1.57 ±0.07 + 0.02 ps 




-> W0uV)K* u 


CDF 


1.42 ±0.09 + 0.02 ps 






CDF 


1.26 ± 0.20 ± 0.02 ps 



Lifetime Ratios (CDF [21], Run II prelim., compared with 
world average (ave) [18] an d HQE predictions [5]): 

■ i |t(B*)/t(B°)=1.11 ±0.09 



— ^ Run I: 1.091 ±0.050 
t- Ave: 1.073 ±0.014 

HQE: 1.06 ±0.02 

1 t(B s )/t(B°)=0.89 ± 0.15 

— i Run I: 0.899 ± 0.072 

Ave: 0.997 ± 0.038 

4; , | HQE: 1.00 ± 0.01 
0.7 0.8 0.9 1 1.1 1.2 



Table 4. Lifetimes from B -> J/if/X at Tevatron Run II. CDF data 
correspond to an integrated luminosity of 70 pb -1 . 



the momentum of the B, which is needed to calculate 
the proper lifetime from r(B) = L xy M(B)/(cp t (B)) is un- 
known. It can however be related to the J/i// momentum 
via 

p t (B) = F(p t (J/if/))-p t (J/if/) 

where F (p t (J /if/)) is the mean ratio p t (B) / p t (J /if/), and 
the uncertainty on p t (B) depends on the the spread of that 
ratio for different momenta. Both the mean ratio and its 
variance are obtained from Monte Carlo. The results of 
such a Monte Carlo study at D0 are shown in Fig. 10. The 
measured average B-hadron lifetimes are 

• D0 (March 03): r B = 1.561 ± 0.024 ± 0.074 ps 

• CDF (July 02): r B = 1.526 ± 0.034 ± 0.035 ps 

which is consistent with the world average of 
r B = 1.573 ± 0.007 ps [18]. 

2.3.5 Semileptonic Decays With £+ IP Trigger 

CDF is also using B — > T>£v X and A b — > K c £v decays 
from the lepton plus displaced track trigger for lifetime 
measurements. The missing v momentum is accounted for 
using the same Monte Carlo-based method as in the inclu- 
sive B lifetime study discussed above. The main challenge 
is to correct the lifetime bias due to the impact parameter 
cuts in the trigger. The acceptance as a function of lifetime 
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Sidebands 
left sideband: 



Mass (L xy > 400 fim) 



MC Information | 



Mean=5.243 ± 0.008 (GeV) 
a = 72 ± 9 (MeV) 
B ± Signal 96 Events 



D0 Run II Preliminary 




5.2 5.4 5.6 5.8 6 
Mass (GeV. c 2 ) 



(Mass resolution has improved by factor 
~ 2 since this plot was produced) 



Signal Region: 



B* Lifetime 




Figure 9. Life time fit at D0 

shown in Fig. 1 1 is found by a detailed Monte Carlo study. 
A fit to the B s — > D sy uv life time distribution is shown in 
Fig. 12. The statistical precision achieved with the cur- 
rent data sample is ct Tb+ = ~ 0.05 ps, <x T = ~ 0.06 ps, 

cr TBs = ~ O.lOps, cr TAb = ~ 0.13 ps. The full results will be 
published as soon as the systematic errors are fully under- 
stood. 

2.3.6 Lifetimes at the Tevatron - Summary & 
Prospects 

The Tevatron is going to provide high statistics samples of 
all B-hadrons, including B s , B c , At,. Preliminary Run II re- 
sults from fully reconstructed hadronic decays are already 
approaching Run I precision, higher statistical precision 
is expected soon from the lepton+displaced track sample. 
The Run Ha projection (MC studies from Dec-01 [3]) for 
the life time ratios are 

• cr(r Bs /T B o), o-(r Ab /r B o) < 1% 

which will provide a real test of theory for the B s and, 
pending improved theoretical calculations, for the Ab life- 
time. 






49.35 / 33 


A 


0.2301 ±0.0194 


B 


-04067 ±0.03085 


C 


0.7S41 ±0.003034 



II 




Figure 10. The F-factor (D0) 



0.50 



c 

QD 
o 

ID 




0.3 
ct(B) (cm) 



Figure 11. SVT acceptance as a function of cr B (CDF) 

Sideband Signal 
CDF Run II Preliminary 70pb 1 B°s-> uvD s 

Signal distribu] 




pseudo-proper decay length (mm) pseudo-proper decay length (mm) 



Figure 12. Fit to lifetime distribution from B s — > D s /iv 
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3 Lifetime Differences 



3.1 Introduction 



The width difference between long and short lived CP 
eigenstates of the B^ d - B^ d system is predicted to be 



^ ~ (9(10)% 



jf ~ 0(1%) 

^ is large enough to be experimentally accessible, soon. 
The width difference is directly proportional to the mass 
difference, 

ATy = A • Am s . 

where the proportionality constant A is ~ 3 • 10" 3 in the 
Standard Model, but the value suffers from large hadronic 
uncertainties [10]. The mass difference Am s , accessible 
through the oscillation frequency in the system, is itself 
an unknown parameter of great interest, and a major mo- 
tivation for installing the precise vertex detectors at CDF 
and D0 during their upgrades for Run II. It is interesting to 
note that a large value for Am s , corresponding to fast B^ os- 
cillations which are more difficult to measure, corresponds 
to a large value for AI^, which makes it easier to measure. 
ATs and Am s are complementary measurements. Given the 
current limits on Am s , a very small value for AI\ would be 
a hint at new physics. 



• Fit lifetime to purely CP-even B° -> D S D S . With 
certain assumptions B^ — > D^D* is predicted to be 
mostly CP even, so that these decays could be in- 
cluded in the analysis. These assumption would have 
to be tested however, for example with a similar an- 
gular analysis as for the B^ — > J/<^0 case. The result 
for the CP even life time can then be compared to the 
mean lifetime from CP-mixed channels to extract the 
lifetime difference. 

• Fit 2 lifetimes to Bg — > J/*/r0. This can have 3 angu- 
lar momentum states, 2 CP even, 1 CP odd. These 
can be disantangled by an angular analysis. 



D^D^. In small 



• The B.R. Method: Assume that the width difference 
is entirely due to CP-even B^ 
velocity (Shifman-Voloshin) limit [3], [22]: 

BR(B° S -> D^Dh = ATs/Ts 



i + ^Ar s /r s 



3.3 Current Values for AT 

Recent results for the width differences in the B system are 



• ^ < 0.18 (95% CL) (DELPHI, 2002) [23] 

• §r < 031 ( 95% CL ) (combined LEP, CDF, for 
1/T 5 = r(B°), using lifetime method) [18] 

• f± = 0.26+^0 (ALEPH, from B.R. method) [22] 



Theory Status 

Recent theoretical predictions for AT^ are (the different re- 
sults are obtained using different expansions of Next-To- 
Leading Order QCD corrections: 



3.4 Prospects for Ar at the Tevatron 

CDF expects the following precisions on AT S by the end of 
Run Ha, from 2fb _1 . The projections assume ^ = 15% 
and are those given in [3] in December 2001. They refer to 
the statistical error only. 



• ff = (2.6!13 ' 10 " 3 I« 

• ff = (3.0+°;*) • 10" 3 [9] using method in [3] 

For Ar^, recent predictions are: 

• ff = ( 8 - 5 ± 2 - 8 ) % t 5 l usin § method in t 10 ] 

• Ma = (9.0 ± 2.8) % [5] using method in [7] 

3.2 Strategies for Extracting Ar^ 

In principle, one could simply fit two exponentials to the 
lifetime distribution of B° decays to mixed CP states. How- 
ever, since ~ (9(10%) only, this method would require 
very high statistics, therefore extra information is needed to 
seperate the CP eigenstates. Possible strategies include [3]: 



• FromB s — > J/if/cp: 

• B s ->D s D s (noD*): 

• B s -> D^D^: 

(assume decay 100% CP even) 

• B.R. method: 
(model dependent) 



^tat(f^) ~ 5% 

^stat(^) ~ 6% 
OStatiff) ~ 2-5% 



Assuming a similar performance for B s — > J/<^0 at D0 we 
arrive at a total statistical uncertainty at the Tevatron for 
2fb _1 of o- stat (ff) ~ 2%, ignoring the B.R. method. A 
more conservative estimate of <x stat ( y^-) ~ 3% is obtained 

if the assumption that B s -> D^D* decays are 100% CP 
even is dropped and decays involving D* are completely 
ignored. 
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4 Conclusion 

Lifetime ratios 

Since they have started data taking, the B -factories have 
brought the error on <t(tb+ /t b o) down to 1.7%, so that the 
experimental accuracy for this ratio is now better than that 
of the HQE prediction. The agreement between theory 
and experiment is very good. Further improvements on 
o~(tb + /tro) can be expected from B-factories and Tevatron, 

d 

soon. 

Large numbers of B^ and Ab are currently being produced 
at the Tevatron. The uncertainty on the lifetime ratios 
<t(tb s /t b o), o~(r\ h I Tgo ) is expected to be below 1% by the 
end of Run Ha. This will provide a real test of HQE for B s 
for which precise predictions exist, while improved theo- 
retical values are needed for for Ab. 

Lifetime Differences 

AT S and Am s are complementary measurements, and both 
parameters combined are sensitive to New Physics contri- 
butions to Bg mixing. Recent calculations predict ^P- = 
9 + 3% [5]. 

From data we get the following limit on the lifetime differ- 
ence in the B s system: ^ < 0.31 (95% CL) [18]. First 

steps have been taken towards a y^- measurement at the 
Tevatron, where 55 ± 9 B s — > J/if/cp events have been re- 
constructed at CDF, and an average B s lifetime has been 
extracted from that decay. By the end of Run Ha a mea- 
surement of y± with a statistical uncertainty of ~ 2% is 
expected. 
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